Regularization of Discrete Ill-Posed Problems

The paper concerns conditioning aspects of finite-dimensional problems arising when the Tikhonov regularization is applied to discrete ill-posed problems. A relation between the regularization parameter and the sensitivity of the regularized solution is investigated. The main conclusion is that the condition number can be decreased only to the square root of that for the nonregularized problem. The convergence of solutions of regularized discrete problems to the exact generalized solution is analyzed just in the case when the regularization corresponds to the minimal condition number. The convergence theorem is proved under the assumption of the suitable relation between the discretization level and the data error. As an example the method of truncated singular value decomposition with regularization is considered. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical differentiation from a viewpoint of regularization theory

In this paper, we discuss the classical ill-posed problem of numerical differentiation, assuming that the smoothness of the function to be differentiated is unknown. Using recent results on adaptive regularization of general ill-posed problems, we propose new rules for the choice of the stepsize in the finite-difference methods, and for the regularization parameter choice in numerical differentiation regularized by the iterated Tikhonov method. These methods are shown to be effective for the differentiation of noisy functions, and the order-optimal convergence results for them are proved.

     

Parameter choice by discrepancy principles for ill-posed problems leading to optimal convergence rates

Schock (Ref. 1) considered a general a posteriori parameter choice strategy for the Tikhonov regularization of the ill-posed operator equationTx=y which provides nearly the optimal rate of convergence if the minimal-norm least-squares solution belongs to the range of the operator (T ^* T) ^v , o<v1. Recently, Nair (Ref. 2) improved the result of Schock and also provided the optimal rate ifv=1. In this note, we further improve the result and show in particular that the optimal rate can be achieved for 1/2v1.The final version of this work was written while M. T. Nair was a Visiting Fellow at the Centre for Mathematics and Its Applications, Australian National University, Canberra, Australia. The work of S. George was supported by a Senior Research Fellowship from CSIR, India.     

2D Burgers equation with large Reynolds number using POD/DEIM and calibration

Model order reduction of the two‐dimensional Burgers equation is investigated. The mathematical formulation of POD/discrete empirical interpolation method (DEIM)‐reduced order model (ROM) is derived based on the Galerkin projection and DEIM from the existing high fidelity‐implicit finite‐difference full model. For validation, we numerically compared the POD ROM, POD/DEIM, and the full model in two cases of Re = 100 and Re = 1000, respectively. We found that the POD/DEIM ROM leads to a speed‐up of CPU time by a factor of O(10). The computational stability of POD/DEIM ROM is maintained by means of a careful selection of POD modes and the DEIM interpolation points. The solution of POD/DEIM in the case of Re = 1000 has an accuracy with error O(10^?3) versus O(10^?4) in the case of Re = 100 when compared with the high fidelity model. For this turbulent flow, a closure model consisting of a Tikhonov regularization is carried out in order to recover the missing information and is developed to account for the small‐scale dissipation effect of the truncated POD modes. It is shown that the computational results of this calibrated ROM exhibit considerable agreement with the high fidelity model, which implies the efficiency of the closure model used. Copyright © 2016 John Wiley & Sons, Ltd.     

Delayed stability switches in singularly perturbed predator–prey models

In this paper we provide an elementary proof of the existence of canard solutions for a class of singularly perturbed planar systems in which there occurs a transcritical bifurcation of the quasi steady states. The proof uses the one-dimensional result proved by V.F. Butuzov, N.N. Nefedov and K.R. Schneider, and an appropriate monotonicity assumption on the vector field. The result is applied to identify all possible predator–prey models with quadratic vector fields allowing for the existence of canard solutions.     

Orientational deconvolution of two-dimensional static disorder by a Tikhonov regularized method for 2H solid state NMR of nano-tubular-oriented structures

A general trend of supramolecular chemistry is the building of mesoscopic-oriented nanotubes. When the typical radius of the channel almost fits the mean supramolecular cross-sectional radius, the guest supramolecule is ultraconfined in the channel. In that case, only rotational disorder around the channel axis is allowed. For sufficiently low temperatures, we expect this 2D disorder to be static on the NMR timescale, and it should reflect the local symmetry of the channels. In this article, we show that experiments performed with deuterium single crystal solid-state NMR as a function of rotation around the channel axis may lead to important information concerning the static orientational disorder of the confined supramolecules. Using an inversion method with Tikhonov regularization and a positivity constraint, that also takes into account of the nonideal pulse sequence response, the orientational probability density can be obtained from 1D experiments performed at different angles, even for significant disorder. As a first step, the method is validated on four different theoretical distributions. It is then applied to ²H NMR single crystal experiments performed on an archetype of parallel channels intergrowth compounds: selectively deuterated 1,10-decanedicarboxylic diacids in hydrogenated urea channels. The diacids form one-dimensional infinite hydrogen-bonded chains ultraconfined in urea linear channels.     

An FEM approximation for a fourth-order variational inequality of second kind

A fourth-order variational inequality of the second kind arising in a plate frictional bending problem is considered. By using regularization method, the original problem can be formulated as a differentiable variational equation, and the corresponding discrete FEM variational equation is presented afterwards. Abstract error estimates and error estimates of the approximation are derived in terms of energy norm and L^2-norm.     

A new zero-finder for Tikhonov regularization

Tikhonov regularization with the regularization parameter determined by the discrepancy principle requires the computation of a zero of a rational function. We describe a cubically convergent zero-finder for this purpose. AMS subject classification (2000)  65F22, 65H05, 65R32